Dynamic control of quantum geometric heat flux in nonequilibrium spin-boson model 
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We study the quantum geometric heat flux in the nonequilibrium spin-boson model. By adopting the 
noninteracting-blip approximation, we show that there exists the geometric heat flux only when the two level 
system is nondegemate. Moreover, the pumping, no pumping, and dynamic control of geometric heat flux 
are discussed in detail, compared to the results with Redfield-weak coupling approximation. In particular, two 
system-bath couplings modulation induced geometric energy tansfer is identified, which is exclusive to quantum 
transport in the strong system-bath coupling regime. 
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Introduction. Smart energy control in low dimensional 
nanoscale systems is of both theoretical and practical impor- 
tance, rendering the emergence of Phononics [1]. As is well 
known, according to the second law of thermodynamics, en- 
ergy cannot transfer from a cold reservoir to a hot side without 
external modulation. In order to obtain a more flexible con- 
trol of thermal energy at nanoscale, there is a great demand in 
designing intriguing phononic devices, which can utilize tem- 
poral modulations to achieve dynamic control, such as in heat 
pump, motor and engine. 

So far, many proposals have been provided to control the 
energy flow between a cold reservoir and a hot one [2-9]. 
As a result of these investigations [3, 6], there is no way to 
pump the energy from the cold side to the hot side in classical- 
oscillator systems by force driving, though heat pump can 
be implemented in classical spin chains [3]. An interest- 
ing design to realize the heat transfer is using the adiabatic 
modulation [10, 11]. In particular, Ref. [10] proposed to 
pump the energy between two reservoirs with equal temporal- 
averaged-temperature by utilizing the geometrc phase induced 
heat transfer. Very recently, similar geometric heat flux in 
classical interacting systems is also studied [11]. 

There are many approximation methods in studying the 
heat transfer in spin-boson model [10, 12-17]. Among 
them, the Redfield-weak coupling approximation is often used 
since this method is effective. However, recent studies re- 
veal that there are also limitations to the Redfield-weak cou- 
pling scheme [15-19]. Besides this scheme, other methods 
such as the multilayer multiconfiguration Hartree [15] and the 
noninteracting-blip approximation (NIBA) [20-22] are also 
applied to study the heat flow under temperature bias. Interest- 
ingly, these methods can produce results different from that of 
the Redfield approximation method for the same Hamiltonian. 
In particular, as shown in Refs. [15, 17], the heat current is not 
linearly dependent on the system-bath interaction strength as 
given by the Redfield-weak coupling scheme. The maximal 
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heat current exists at the intermediate system-reservoir cou- 
pling strength[15, 17], and then the heat current descends as 
increasing the system-reservoir coupling strength. So far, dif- 
ferent physical results of Redfield scheme and NIBA are not 
found in the weak system-bath interaction regime. 

There are mainly two assumptions inherent in the Redfield- 
weak coupling scheme. One is the resonant tunneling be- 
tween the system and the reservoirs, and only the resonant 
frequency in the reservoir contributes to the dynamics. The 
other is that each individual reservoir interacts with the system 
separatively, or say, additively. In this paper, we shall not use 
these two assumptions, but adopt the NIBA scheme to study 
the geometric phase induced energy transfer The effect of 
non-resonant tunneling between the system and the reservoirs 
as well as their collective non-additive interacting are both 
considered. We shall first review the analytical expressions 
of the investigated spin-boson model in NIBA, and the geo- 
metric phase induced energy transfer through the generating 
function approach. We then investigate in detail the behaviors 
of geometric heat flux under various modulation protocols. 

Our contributions are two-folds. First, we find different re- 
sults of the two methods (Redfield and NIBA) with whatever 
system-bath interaction strength through calculating the geo- 
metric phase induced heat flux in unbiased (degenerate) case. 
This renders the geometric heat flux could become an effective 
tool to compare the difference or even judge the correctness of 
various approximation schemes. Second, geometric heat flux 
itself is a very important physical problem. A more thorough 
investigation to the problem with different methods is mean- 
ingful. We find that the temperatures-modulation-induced ge- 
ometric heat flux are absent under the unbiased case, while 
linearly dependent on the system-bath coupling strength under 
the biased (nondegenerate) case. Also, additional to reversing 
the modulation protocol, the direction of geometric heat flux 
can be reversed by adjusting the system's parameters. More- 
over, two system-bath couplings modulation induced geomet- 
ric heat flux is identified. So far, this nontrivial observation 
is exclusive for quantum transport in the strong system-bath 
coupling regime. 

Model. The non-equilibrium spin-boson (NESB) model, 
[see the left panel of Fig. 1], consisting of a two-level sys- 
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FIG. 1: (color online). Schematic illustrations of the nonequilib- 
rium spin-boson model for energy transfer. The left panel describes 
the Hamiltonian (1) while the right one depicts the Hamiltonian (2). 
They are equivalent, with only a rotation difference of the basis. 



tern in contact with two bosonic reservoirs with temperatures 
(z/ = L, R), is described by the Hamiltonian: 



H = — (T7 
2 



(1) 

where is the energy gap of the two levels; A denotes the 
tunneling between them; = — |0)(0| and = 

|0)(1| + |1)(0| are the Pauli matrices expressed in the two 
level basis; foj ^{bj,v) denotes the creation (annihilation) op- 
erator of the jth harmonic mode in the v bosonic bath, with 
Xj,,y the system-bath coupling strength. We can also rotate the 
Hamiltonian (1) firstly around the y axis by an angle 7r/2, 
then around the new x axis by an angle tt. In this way, 
cTx "> CTzjO-z <^x, and the alternative equivalent Hamil- 
tonian reads, [see the right panel of Fig. 1]: 



A 



(2) 

Throughout the paper, we will work in the basis of Hamil- 
tonian (1) in the NIBA scheme. Before proceeding to the en- 
ergy transport, it is useful to transform the NESB Hamiltonian 
(1) by using the canonical transformation [23] (also called 
Lang-Firsov or polaron transformation): Hx = HU, U = 
exp[ia,n/2], n = 2iJ2^^j ^(^> - ^j,'^)- After neglect- 
ing an irrelevant constant u ^^'j i^/'^i' '■^^ transformed 
Hamiltonian is expressed as: 

A 



IT '^0 



with CT+ = {(Tx + iay)/2 = |1) (0| and cr_ = {<Jx — iay)/2 = 
|0)(1|. Ht clearly shows that the energy transfer is accom- 
plished by the excitation from the lower level to the upper one 
with absorbing energy from the baths, and the relaxation from 
the upper level to the lower one with releasing energy to the 
baths. Thereafter, as a result of the NIBA method [19-22] and 
with the Markov assumption, the population dynamics of the 
two-level system becomes [17]: 
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dt 
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where, Po/i(0 = (1 T denotes the population at 

the lower (upper) level. The transition rates stand for the co- 



operative process between the system and two reservoirs: 

(A/2)2 1-^ 



2tt 



C^{uj)duj, 



(5) 



C±H = Ci(±£oT^)C,j(±c^), 



with Cu{uj) = /^e'"*~'3.(t)rft denoting the proba- 
bility density of the bath v to absorb the energy lo 
(equivalently, to release the energy —uj). Employ- 
ing the Gaussian statistics of the "momentum" opera- 
tor we have Q^{t) = ([17(0) - ^2(01^^(0)) = 
/o°° [coth(2f^) (1 - cos(wt))+i sin(wt)]dw [22], with 
Jj/(w) = 47r A|^(5(w — (jJj^u) being the spectral density of 
the bosonic bath v. In contrast to the Redfield-weak coupling 
scheme [10, 12], the rate expressions (5) distinctly exhibit the 
non-resonant energy tunneling processes, conjoining the two 
baths nonadditively: K{eQ) means when the central system 
loses energy Eq by relaxing from the upper level to the lower 
one, the R bath will absorbs uj and the L bath gains the rest if 
Eq > cli or even supplements the shortage if a; > Eq, K{—£q) 
depicts the similar energy transfer process for the central sys- 
tem exciting from the lower level to the upper one. 

Generating function and geometric heat flux. Following 
the Full Counting Statistics [24, 25], we next construct the 
cumulant generating function (CGF) of the NESB model to 
count the phonon energy transfer through the right system- 
bath coupling [10, 11, 18, 26-28]. Denote pt{n,uj) as the 
joint probability that a total energy of uj has been transferred 
to the right bath during time interval [0,t], with the two-level 
system populates at state \n) (n — 0, 1) at time t, we then 
introduce the characteristic function of that joint probabil- 
ity |z(x,/)) ^ U°°^pt{Q,w)e^-^dL,,!^^p,{l,uj)e^-^du:Y 
with the energy counting field x- Following [17], this charac- 
teristic function satisfies the following dynamics: 



z{x.t))=-'H{x)\z{x,t)), 



dt 

with ■H{x) = 



K{-eo) -K+ix) 
-K-{x) K{eo) 



(6) 



where i^=^(x) = C±(w)e±*">fda;. When the 

counting field x = 0, K^{0) = K{zLeo) and in turn Eq. 
(6) reduces to Eq. (4). Thus, the characteristic function of the 
heat transfer is Z{x,t) = /^^ (pt(0,w) + pt{l,u;))e''^^duj 

and the CGF is G{x) = linij-i-oo 7 ln[^(x, 01' which gener- 
ates the n-order cumulant of heat transfer fluctuations through 
limt^ooiiQ")) /t = d"^g{x)\x=o- The mean value of the 
heat flux is just the first order: J — di^Q (x) |^=o 

Behaviors in the long time limit are of our central inter- 
est. They are governed by the ground state of the operator 
'H(x), of which the eigenvalue i?o(x) possesses the the small- 
est real part. For time-independent H(x), linit-s-oo Z{x, t) ^ 
e-Eo(x)t and in turn linit^oo((Q"))/t = -dll^Eo{x)\x=o- 
However, for time-dependent H{x,t), where the system 
parameters A{t), eo{t), the bath temperature T^{t) or the 
system-bath coupling Xj^ui^) could be subject to periodic 
modulations, the adiabatic geometric phase effect has been 
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unravelled to play an important role in the dynamic control of 
energy transfer [10, 11]. In this case, although at every instant 
the dynamics (6) is preserved, there exist two contributions 
intheCGF: limt^oo ^(x, ^ e*^ = gt(e<i„„+egeo™)^ One 
is the dynamic part Gdyn, and the other is the geometric part 

Ggecn [10, 11,28]: 



Qdyn — 



1 



dtEoix,t), 
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dt{(l>o\dt\ipo), 



with Tp the modulating period [29] and {{(j)o\) the bi- 
orthonormal right (left) eigenvector corresponding to the 
ground state of 'H{x, t) [30]. In the case of two parameters 
being modulated, say ui{t),U2{t), the calculation of Ggeom 
can be done using Stokes theorem [31], 



duidu2Buii 



(8) 



and the Berry curvature [31, 32] is: 

Bmu2 = (c'„i(/)o|9„2?Ao> - (9„2 0o|c^«iV'o) 

_ (0o|gm'H|^l)('/)l|9M2"^|V'o) " (ui O U2) 



(9) 



with El the eigenvalue of the excited state and \tpi) {{(j>i\) 
the corresponding bi-orthonormal right (left) eigenvector. Ac- 
cording to these formulas, we have the dynamic heat flux Jdyn 
and the geometric heat flux Jgeom, 



J dyn 



Jo. 



1 

Tp JQ 



Tp 



d{ix) 
duidu2 



d{ix) 



(10) 
(11) 



x=o 



Results and discussions. With equations above, we are 
ready to study the consequences of dynamic control of the 
NESB model in a more exact way, i.e., without Redfi eld- weak 
coupling approximation. 

i) Unbiased case, eo = 0. In this degenerate case, we have 
is:(-eo) = K{eo),K+{x) = K-{x), and the matrix. 



n{x.t) = 



K{eo) -K+ix) 
-K+ix) K{eo) 



(12) 



Thus, the corresponding eigenvalues and eigenvectors turn 
into a simple form. 



Eo/i^Kieo)TK+{x), 



I ■00/1 
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'7^ 



±1, 1 



V2 



±1, 1 



Since the eigenvectors are constant, we have that for whatever 



two-parameter modulations, the Berry curvature Bu^ i 



0. 



There is no geometric phase effect and geometric heat flux is 



absent when the two-level energy gap is zero. This result is 
independent of the system-bath coupling strength. 

Note that, in contrast to the absence of geometric heat flux 
here, Ref. [10] treats the same physical system with Eq — 0, 
but in the Redfield-weak coupling scheme and in the basis 
of Hamiltonian (2) (with notation A — > ojq). Nonzero geo- 
metric heat flux is investigated there. The difference proba- 
bly comes from two folds: 1) the effect of the non-resonant 
tunneling and collective nonadditive interacting between the 
system and baths in NIBA scheme, are inherently overlooked 
in the Redfield-weak coupling approximation; 2) the rotation 
difference of the basis. Generally, the different choice of basis 
should not lead to different physical results. But in order to get 
the final population dynamics, the secular approximation are 
applied for both methods, to decouple the diagonal and non- 
diagonal reduced density matrix elements [33]. In this sense, 
after neglecting the off-diagonal dynamics, the two popula- 
tions in the upper/lower level in the rotated Hamiltonian (2) 
[Po/i{t) = (1 T (o'a;(i)))/2] are different from the popula- 
tions in the Hamiltonian (1) [po/i{t) = (1 =F {az{t)))/2]. To 
a final conclusion of the difference, a detailed study of the full 
reduced density matrix dynamics is in demand. 

Any way, although it has been known already that the dy- 
namic heat flux for the same Hamiltonian calculated from the 
two methods are different for intermediate coupling strength 
and essentially the same in the weak coupling limit [17], here 
we for the first time show the difference for whatever cou- 
pling strength with geometric heat flux. Since the geometric 
phase effect contains the high-order heat transfer fluctuations, 
we expect that the geometric heat flux could be an effective 
tool to compare the difference or even judge the correctness 
of various approximation schemes. 

ii) Biased case, Eq ^ 0. An applied magnetic field 
could control this Zeeman splitting. In this nondegenerate 
case, there do exist the geometric heat flux. In the fol- 
lowing, to simplify the calculation we assume Marcus limit 
[22, 34] that works in high temperature T,y > Eq and/or 
strong system-bath coupling regime. The Marcus limit could 
be approached by a short time expansion of Qi, (t) such that 
Q,y{t) — Ti,T^P + iTjyt with the renormalized system-bath 
coupling = / "^^^duj = J2j 4A^^/a;j,^. In this way, we 
have the transition rates [17, 22]: 



Ki±Eo) 



A 

A2 

4 



exp 



FlTl + TrTr 



exp 



(£0 tTlT Tfl)^ 
'4{TlTl+TrTr) 



where A, eo, T^, could be subject to the dynamic control. 
By substituting these rates into the dynamics. (6), we are able 
to investigate the Berry curvature and the geometric heat flux 
Jgeom through Eqs. (9) and (11). 

We first adiabatically modulate the two bath temperatures. 
The external control frequency Hp is chosen to be sufficiently 
small so that the adiabatic condition is respected [29]. The 
protocol is chosen as TlI^) = 150+90 cos(l]pt) (K),Tfl(t) = 
150 + 90sin(ript) (K), so that there is no temperature-bias- 
induced flux on average (Jdyn — 0), but the geometric heat 
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FIG. 2: (color online). 
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Integrated geometric heat flux per period 
Qgeom = JgaomTp as a linear function of the system-bath coupling 
under the two-bath-temperature modulation. We set the symmetric 
coupling r_L = Vr — r and tunnelling energy is A = 5.22 meV. 
When eo — eo, we observe the same lines. 



flux emerges. Figure 2 shows that the integrated geometric 



heat flux per period Qg 



Jgeomlp IS Unearly depen- 



dent on the system-bath coupHng F. When F approaches to 
zero, Qgeom does not vanish but persists at some finite values. 
Also, we find that increasing the gap Eq can increase Qgeom 
and when Eq — > — £o, we observe the same lines. Generally, 
reversing the modulation cycle can reverse the direction of ge- 
ometric heat flux. Moreover, as indicated in Fig. 2, when the 
coupling strength F exceeds some threshold values the geo- 
metric heat flux will also reverse its direction. This may offer 
a useful mean in dynamic control of the heat flux induced by 
adiabatic periodic modulation. 

Second, we manipulate the bath temperature and the two- 
level energy gap Eq- As shown in Fig. 3(a), Qgeom is lin- 
early dependent on the system-coupling strength. When F ap- 
proaches to zero, Qgeom vanishes. Raising the average level- 
gap Eg is able to increase the magnitude of Qgeom- If we re- 
verse £o the geometric heat flux reverses its direction 
so that when Eq = 0, Qgeom is absent. Although in this mod- 
ulation protocol the dynamic heat flux is nonzero, it decays 
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FIG. 3: (color online), (a) Qgeom as a linear function of V un- 
der modulations of one bath temperature and the gap eo (i) • We set 
TL{t) = 150 + 90 cos(J7pt) (K),eo(t) = £o + 0.78sin(ript) (meV). 
Tr = 150 K, A = 5.22 meV. (b) Integrated dynamic heat flux per 
period Qdyn = JdynTp under the same conditions for comparsion. 
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FIG. 4: (color online). Emergence of geometric phase effect and 
geometric heat flux for modulating two system-bath couplings. Al- 
though Qgeom is absent at Tl ~ Tr, the nonzero geometric skew- 
ness (d'i^Ggeom\x=o) shows the existence of the geometric phase ef- 
fect, manifesting itself as the high-order heat transfer fluctuations. 
The control protocol is Tl = 130.5 + 104.4 cos(Slpt) (meV) 
and Tr = 130.5 -I- 104.4 sm{flpt) (meV). Other parameters are 
Tl = 150 K, a;o = 2.61 meV and A = 5.22 meV. 



very fast as the system-bath coupling F increases, [see Fig. 
3(b)]. In this way, the geometric heat flux will dominate the 
energy transport in the strong system-bath coupling regime. 

Besides these two control protocols discussed above, we 
finally consider the special case of modulating two system- 
bath couplings ^^{t) and Ffl(t). Figure 4 shows the emer- 
gence of geometric heat flux when ^ Tjj. Although at 
Tl = Tff, the geometric heat flux vanishes, the geometric 
phase effect still exists, manifesting itself as the high-order 
heat transfer fluctuation, e.g., the nonzero geometric skewness 
dfy^Ggeom\x=Q 7^ 0. Thcsc observations are distinct from the 
previous results either in quantum weak-coupling regime [10] 
or classical regime [11], wherein under couplings-modulation, 
the BeiTy curvatures are always zero (no matter what other pa- 
rameter settings), so that the geometric phase effect and geo- 
metric heat flux are always absent. It is the strong system- 
bath coupling we consider here that makes the couplings- 
modulation-induced geometric heat pump nontrivial. 

In summary, we have studied the geometric phase induced 
heat flux extensively in the spin-boson system under the adia- 
batic periodical modulation without the Redfield approxima- 
tion. Using NIBA approach, we have found that there ex- 
ists the geometric heat flux only when the two level system's 
energy gap is not zero. This result shows that NIBA and 
the Redfield approximation can give different results for the 
same physical system even when the coupling is weak. More- 
over, the pumping, no pumping, and dynamic control of ge- 
ometric heat flux have been discussed in detail. In particu- 
lar, two system-bath couplings modulation induced geometric 
heat flux has been identified. So far, this nontrivial observa- 
tion is exclusively for quantum transport in the strong system- 
bath coupling regime. 
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